Field-induced phase transitions of repulsive spin-1 bosons in optical lattices 
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We study the phase diagram of repulsively interacting spin-1 bosons in optical lattices at unit 
filling, showing that an externally induced quadratic Zeeman effect may lead to a rich physics 
characterized by various phases and phase transitions. We find that the main properties of the 
system may be described by an effective field model, which provides the precise location of the 
phase boundaries for any dimension, being in excellent agreement with our numerical calculations 
for one-dimensional (ID) systems. Our work provides a quantitative guide for the experimental 
analysis of various types of field-induced quantum phase transitions in spin-1 lattice bosons. These 
transitions, which are precluded in spin-1/2 systems, may be realized using an externally modified 
quadratic Zeeman coupling, similar to recent experiments with spinor condensates in the continuum. 
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Ultracold atoms in optical lattices constitute a 
highly controllable scenario for the analysis of strongly- 
correlated systems, as highlighted by the realization of 
bosonic and fermionic Mott-insulators (Mis) [H-Q ■ Inter- 
estingly, on-going experiments 0, HI are approaching the 
regime at which magnetic properties, including the long- 
pursued Neel antiferromagnet in spin-1/2 fermions, could 
be revealed. Optically trapped spinor gases, formed 
by atoms with various Zeeman substates, are particu- 
larly interesting in this sense. The internal degrees of 
freedom result in an exceedingly rich physics, mostly 
studied in the context of spinor Bose-Einstein conden- 
sates (BECs) [6|-l9|]. Spinor gases in lattices are particu- 
larly exciting, since they provide unique possibilities for 
the analysis of quantum magnetism. 

Spin-1 gases are the simplest spinor system beyond the 
twcHComponent one. Depending on interparticle interac- 
tions [SIlZI (given by the s-wave scattering lengths ao,2 for 
collisions with total spin and 2), spin-1 BECs present a 
ferromagnetic (FM) ground state (for ao > a-i as in 87 Rb 
F = 1 19(]) or an antiferromagnetic (AFM), also called po- 
lar, one (for a 2 > a , as in 23 Na [8j). Spin-1 lattice bosons 
have also attracted a strong interest, especially the AFM 
case, for which a wealth of quantum phases have been 
predicted [Tol4T8|] . For AFM interactions, in 2D and 3D 
the MI states at odd filling are nematic 13, 17, Eij ]. 
whereas in ID quantum fluctuations lead to a sp onta- 

noeusly dimerized ground state dEEl EE BSE The 

case ao = ai exhibits an enlarged SU (3) symmetry with 
a highly degenerate ground state [23[ . 

Most spin-1 species are naturally close to this SU(3) 
point (i.e. a w 02), where small external perturba- 
tions, as Zeeman shifts, may have a large effect, re- 
ducing the system symmetry, and thus favoring differ- 
ent phases. Since interactions preserve the magnetiza- 
tion, Ai, the linear Zeeman effect (LZE) may be effec- 
tively gauged out (although the phase diagram depends 
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FIG. 1: Mott phases of spin-1 lattice bosons at unit filling, 
as a function of 9 = arctan(J2/Ji) — tt and the QZE D. Thick 
solid lines correspond to first order phase transitions for any d. 
Grey region is the dimerized phase (only in ID). The symbols 
represent ID numerical data (see text). 

on M. [13, EH)- On the contrary, the quadratic Zee- 
man effect (QZE) plays a crucial role in spinor gases. 
In spite of its importance, the role of the QZE in the 
quantum phases of spin-1 lattice bosons remains to a 
large extent unexplored, with the sole exception of the 
recent 3D mean-field analysis of Rcf. 



18[, where it was 



shown that for finite M the QZE may lead to nematic- 
to-ferromagnetic (or partially magnetic) transitions. 

This Letter discusses, for the first time to our knowl- 
edge, the complete phase diagram (Fig. [IJ for MI phases 
(with unit filling) of spin-1 bosons in the presence of QZE, 
for the experimentally relevant case of a balanced mix- 
ture, i.e. with M. = 0. Combining an effective field the- 
ory (for any dimension), ID DMRG calculations, and ex- 
act Lanczos diagonalization, we obtain the phase bound- 
aries, characterizing the phase transitions. We note that 
the QZE may be controlled by means of microwave and 
optical techniques [2^, [26[. Hence, as recently demon- 
strated for spinor BECs in the continuum [24| . our re- 
sults show that a controlled quenching of the QZE may 
permit the observation of field-induced phase transitions 
in spin-1 lattice bosons, which are precluded by simple 
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use of the LZE due to conservation of M., and thus are 
absent in spin- 1/2 systems . In addition, optical Feshbach 
resonances [2?], [28| permit the modification of the ratio 
0,2/0,0, so that the full phase diagram discussed below 
may be explored with state of the art techniques. 

We consider repulsively-interacting ultracold spin-1 
bosons in a d-dimcnsional hypercubic lattice, prepared 
in a balanced mixture {Ai = 0). The inter-particle in- 
teractions are characterized by the coupling constants 
.90,2 = 47r?i 2 ao l 2/w a (where m a is the atomic mass). At 
integer filling the system is in the MI regime if the (pos- 
itive) interaction constants 170.2 3> t, where t is the hop- 
ping amplitude between neighboring sites. In second or- 
der of perturbation theory in f, the low-energy physics is 
given by super-exchange processes, being described by an 
effective bilinear-biquadratic spin Hamiltonian 



jea by a 
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Hr=- 



J\Si ■ Sj + J2\Si ■ Sj 



(1) 



where Si are spin-1 operators at lattice site i, the sum 
runs over nearest neighbors, J\ = 2t 2 /g2 and J2 — 
2t 2 /3g 2 + 'it 2 /3g (both are positive). Note that the FM 
case (ao > a 2 ) corresponds to J\ > J2, whereas the AFM 
case (<Z2 > ao) results in J2 > J\. As mentioned above, 
typically ao « a 2 , which corresponds to the vicinity of 
the ferromagnetic SU(3) point (Ji = J2). 

As aforementioned, the QZE, characterized by the ex- 
ternally controllable constant q, plays a crucial role in 
the physics of the system. The QZE (which resembles 
the so-called single-ion term in other condensed-matter 
scenarios) leads to a modified Hamiltonian of the form 



h = h i -djY,{si)\ 



(2) 



where D = q/J, with J = \/ J 2 + Jf . In the follow- 
ing, we employ J as the energy unit, setting J = 1, and 
introduce the standard parametrization of the exchange 
constants J\ = — Jcos(9), J 2 = —Jsin(0), where the 
angle 9 takes values in the interval (— it + arctan |, — f ) 
as the ratio 02/00 varies from to +00 [13(. 

For 9 < -37r/4 ( J x > J 2 ) if D > the ground state is a 
fully polarized ferromagnet in either m = ±1 (Ising-FM) , 
and since M. = phase separation into ferromagnetic 
to = ±1 domains is expected. For D < the ground 
state for small values of \D\ is an XY- ferromagnet (XY- 
FM), i.e. the system fulfills (Sf) = 0, but presents a non- 
zero transversal magnetization. This phase is ordered in 
dimensions d > 2, exhibiting in ID a quasi-long-range or- 
der with leading power-law decay of xy spin correlations. 
For larger \D\ (keeping D < 0) there is a phase transi- 
tion between the XY-FM and the so-called large- D phase 
(also called Ising-nematic), in which all atoms are in the 
to = Zeeman substate, and hence all spin correlations 
decay exponentially. The field-induced phase transition 
between large- D and XY-FM is discussed in detail below. 



the dominant correlations 
ll,[2l|. A XY- 



For 9 > -3tt/4 (Ji < J 2 , 
are of spin-nematic (quadrupolar) type 
nematic phase occurs for D > 0, characterized for d > 
2 by ((S+) 2 ) ^ and (S) = 0, and in ID by power- 
law correlations of the quadrupolar order parameter and 
exponentially decaying in-plane spin correlations. On the 
contrary, for D < the large-D phase is favoured. In 
ID, for D = a dimer nematic phase is expected 



m llS 12J, l22|- We show below that in ID the QZE 



induces a transition between the dimer-nematic phase 
and the XY-nematic (or large-D) one. 

To study the phase diagram near the SU(3) point, we 
develop a low-energy effective field theory [291 ] based on 
spin-1 coherent states \ip) = J2a=x y z (. u a + iva)\ta) where 
\t a ) are three Cartesian spin-1 states (\t z ) = \m = 0), 
|m = ±1) = T(l/V2)(\tx) ± i\t y )))- The real vectors u, 
v (defined at each lattice site n) satisfy the constraints 
u 2 +v 2 = 1, u-v = 0. The vector u plays the role of 
director vector for the nematic phases discussed below. 
The average spin on a site fulfills M = (V'|S|V') = 2(ux v) 

The FM region 9 < —37r/4 is characterized by M and 
the Euler angles ($, tp, \) which parametrize the orienta- 
tion of the mutually orthogonal vector pair (M, u). The 
Euclidean Lagrangian for the model ^ takes then form: 

L = iJ2M n {drXn + (dr<Pn)cOsK}-Dj2K Z 
n n 

- {jiM n -M n ,+J 2 tr(B£-B n ,)}, (3) 

(nn ! ) 

where B ab = (ip\S a S b \^A is related to the nematic ten- 
sor Q = \ — B — B r I30j| and the magnetization M a = 
tabcB ■ Assuming u and v as smooth fields, we may then 
pass to the continuum, performing a gradient expansion 
in Eq. (J3J. For D < 0, configurations with $ ~ ir/2 and 
X — are favoured. Since (cos 1}, ip) and (M, \) are pairs 
of conjugate variables, and \ become "slaves" and can 

be integrated out (see Ref. [29j for details). The potential 

1/2 

energy is minimized at |M| = Mq = {l — rj 2 } , where 
i] = |Z3|/2Z(sin6' — cos6>), with Z the lattice coordina- 
tion number. We expand the effective action around the 
equilibrium value, |M| = Mo + 5, integrate out S, and 
obtain the effective action for ip, which can be cast in the 
familiar form of the (d+ l)-dimensional XY model 



Axy - (A d - 1 /2g) I d d+1 x(d^) 2 



(4) 



Here A is the ultraviolet lattice cutoff, and the coupling 
constant g, acting as an effective temperature, reads 
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where (S 2 ) = g s C d T: 1 - d f* dk k*- 1 / ^mj + k 2 is the fluc- 
tuation strength, g 2 — (8Zrj 4 )/(\ + 2rf) and m 2 — 
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Z(l — T] 2 )/(X + i] 2 ), respectively, the coupling constant 
and mass of the 8 fluctuations, A = tan (9/(1 — tan#), and 
= (4n) d / 2 T(d/2). 
Model @ describes a phase transition between XY-FM 
and large- D occuring at a non-universal g — g c . For d — 
1 this is a Kosterlitz-Thouless (KT) transition, and the 
XY-FM phase has only a quasi-long-range order. For d > 
2, the phase transition belongs to the (d+ l)-dimensional 
XY universality class, the XY-FM phase is ordered with 
a spontaneously broken U(l) symmetry (ip = tpo), and 
the order parameter (cosipoS x + sin <foS y ) ^ 0. Once 
g c is known, Eq. ([5]) constitutes an implicit equation to 
determine the transition curve D(9). The transition line 
has a universal slope for A — > oo (SU (3) point) given 
by r] = 1 — 0(A -1 / 2 ). Fig. Q] shows the curve for d = I 
obtained after adjusting the single fitting parameter g c 
to the numerical results discussed below. An excellent 
agreement with the numerics is obtained for g c 0.6. 

We have numerically evaluated the large- -D to XY-FM 
boundary in ID by means of DMRG calculations (fol- 
lowing the method of Ref. [3l| for up to 42 sites). We 
found that the most efficient way to locate the phase 
boundary is to study the fidelity susceptibility (FS) [32j ]. 
X(D) = -21im A ^oln{|^(£>M£> + A£>))| 2 }/AD 2 , 
where the quantity under the logarithm is the fidelity, 
i.e. the Hilbert-space distance between the ground states 
at two values of the QZE coupling, D and D + AD. Fig- 
ure [2^, shows the evolution of the peak in x(-D) with 
increasing system size. The finite-size scaling of the peak 
position as a function of the number L of sites follows 
very accurately a 1/L 2 law, confirming its KT character. 
Extrapolating the peak position to L — oo yields the 
curve shown in Fig.[TJ which, as mentioned above, agrees 
perfectly with the effective field-theoretical description 
after fitting the single parameter g c . 

For d > 2, g c may be estimated by neglecting fluctua- 
tions of M and demanding the critical \D\ to match the 
Ising value ZJ\ at J2 = 0, which yields g c — (8Z/5) 1 / 2 . 
For a square lattice the resulting g c 2.53 compares fa- 
vorably with the known Monte-Carlo result g c ps 2.20 
for the classical 3D XY transition on a cubic lattice [33[ . 
The corresponding transition curves obtained from the 
implicit Eq. ([5]) are also shown in Fig. [TJ 

For the AFM region, 9 > — 37r/4, the effective th eory 
can be formulated in terms of the nematic director u [29j j , 
or alternatively via the nematic tensor Q ab = u a Ub—8 a b/3 



•An 



(A d -74 5 „) d d+1 x{(d,Q ab ) 2 + 2A 2 m 2 n Q zz } (6) 



where the coupling g n = J Z/2X vanishes at the SU (3) 
point Q and m 2 = D/ siaff is the QZE-induced mass. 
The u anisotropy is of the easy-plane (easy-axis) type 
for D > (D < 0). For d > 2 there is a long range 
nematic order for any D, with a single transition at D = 0, 
between a gapless XY-nematic phase with ((S + ) 2 ) ^0 at 
D>0 and a gapped Ising-nematic ((S z ) 2 ) = at D<0. 
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FIG. 2: a) Evolution of the fidelity susceptibility x(D) at 
the boundary between XY-FM and large-D phases, with in- 
creasing system size L; b) Chirality r as a function of D, for 
two cuts at the FM (9 = -0.767r) and nematic (9 = -0.547r) 
sides. 

For d — 1 , at D — the AFM phase presents exponen- 
tially decaying nematic correlations (albeit with a very 
large correlation length £ cx e 4 ^^ 9n and a tiny spin gap 
A oc £ _1 close to the SU(3) point [13,[2l|). In this phase 
the lowest excitations have total spin 5 = 2, which can be 
understood by noticing that in model © the magnetiza- 
tion M is a composite field, M a cx e abc Q M {d T Q cd ). Q ab 
can be approximately described as a free field with finite 
mass m = A 1 / 2 , so its correlator corresponding to S = 2 
excitations decays as x~ 1 ^ 2 e~ mx , while the spin correla- 
tor describing 5 = 1 excitations decays much faster, as 
x~ 2 e~ 2mx . However, as aforementioned, this featureless 
disordered-nematic phase [2(| [2l| acquires a very weak 
long-range dimer order all the way to the SU( 3) point, 
due to the condensation of Z2 disclinations 12|, |34[. Al- 
though model ([6]) cannot describe dimerization, it may 
be employed to determine the boundaries of XY- Nematic 
and Ising nematic phases at very small \D\, provided that 
Z2 disclinations do not play role at the corresponding 
phase transitions. When D > 0, this transition is KT, 
and the transition line close to the SU(3) point fulfills 
D = CJ 2 exp{-47r(I - Ji/J 2 )" 1/2 }, where C > is a 
numerical constant. At D < the phase transition is 
Ising-like, with a boundary similar as for D > 0, but a 
negative constant C. As shown below, this analysis pro- 
vides a good insight on the dimer-to-nematic transitions. 

To characterize numerically the boundaries of the 
dimerized phase (where FS remains featureless), we em- 
ployed level spectroscopy analysis [35j]. In a finite chain, 
two dimerized ground states (degenerate in the ther- 
modynamic limit) split in energy, so the lowest excited 
state in the dimerized phase is unique and belongs to 
the A4 = sector. In contrast, both in the large-D 
phase {D < D~) and in the XY-nematic (D > D+), 
the lowest excited states are twofold degenerate, having 
A4 = ±1 and M. = ±2, respectively. Thus, in finite 
chains a level crossing between the lowest excited sin- 
glet and doublet states occurs when changing D. Our 
extrapolated results for Df, obtained by Lanczos diag- 
onalization for periodic systems of up to L = 16 sites, 
are shown in Fig. [T] with A -> D+, □ -» D~. Note that 
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when approaching the SU (3) point our numerics cannot 
recover the exponentially small dimer region, which basi- 
cally reduces to D = line. The finite size extrapolation 
of Z?+ follows 1/L 2 law, confirming its KT nature. 

To determine the universality class of the D = D~ 
transition, we have computed the central charge at D = 
0, 9 = — 0.73-7T 36j. The block entanglement entropy 



for an open ID system of size L, divided into two pieces 
of size I (block) and L — I (environment), behaves as 
5=| log [— sin (77)] + A, where c is the central charge 
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Setting / = f 



and A is a non-universal constant 

following Ref. (39|, and using DMRG to evaluate S for 
several L values, we obtain c ~ 1.5. The D = KT 
line has c = 1; subtracting its contribution, we get c = i 
for the D = D~ line, confirming its Ising nature. 

Finally, we discuss the behavior of the chirality r = 
T Eifo+i + <-i - 2 <o) = \ {3(Sf ) 2 - 2}, which 
can be easily monitored in Stern- Gerlach-like time-of- 
fiight experiments. Figure [2}d shows r as a function of 
D. At the FM side, r is discontinuous at D = indi- 
cating the first-order character of the transition, which is 
clear since for the Ising-FM phase (D > 0) r = 1, while 
in the XY-FM phase (D < 0) for D -> -0 the ground 
state energy is minimized at (M, 1?, %) = (1,§,0), and 
thus for ((S z ) 2 } -> I and r -> -±. At the XY-FM to 
large-D transition, r practically saturates to —2 for a 
value of D close to that obtained from the FS analysis. 
On the AFM side the limit D — > is non-singular and 
thus r = there. The nematic-to-dimer transitions do 
not present any pronounced feature of r. These transi- 
tions could be revealed experimentally by Faraday rota- 
tion techniques [40( or those recently explored in Ref. 41 1 . 

In summary, we have obtained the complete phase dia- 
gram (for any dimension) for spin-1 lattice bosons in the 
MI phase (at unit filling) in the presence of quadratic 
Zeeman coupling. Our results provide hence a quantita- 
tive guide for the analysis of field-induced quantum phase 
transitions in lattice bosons, which, similar to recent ex- 
periments with spinor BECs in the continuum [jij ]. may 
be realized modifying the QZE by means of microwave 
dressing. Starting in the large-D phase, and dynamically 
modiying the QZE across the transitions discussed in this 
paper, should result in the FM regime in the appearance 
of XY-FM domains, similar to those observed in spin-1 
BECs [24j , whereas quenches in the AFM regime should 
lead to nematic domains with different ((S x , y ) 2 ) but ho- 
mogeneous (S) = 0. We stress that such field-induced 
transitions are precluded for spin-1/2, constituting an 
interesting novel feature of lattice spinor gases. 
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